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ABSTRACT 

Young Stellar Objects (YSOs) and in particular protostars are known to show a variety of 
high-energy processes. Observations in the X-ray and centimetric radio wavelength ranges are 
thought to constrain some of these processes, e.g., coronal-type magnetic activity. There is a 
well-known empirical correlation of radio and X-ray luminosities in active stars, the so-called 
Giidel-Benz relation. Previous evidence whether YSOs are compatible with this relation remains 
inconclusive for the earliest evolutionary stages. The main difficulty is that due to the extreme 
variability of these sources, simultaneous observations are essential. Until now, only few YSOs 
and only a handful of protostars have been observed simultaneously in the X-ray and radio 
range. To expand the sample, we have obtained such observations of two young clusters rich in 
protostars, NGC 1333 and IC 348. While the absolute sensitivity is lower for these regions than 
for more nearby clusters like CrA, we find that even in deep continuum observations carried out 
with the NRAO Very Large Array, the radio detection fraction for protostars in these clusters is 
much lower than the X-ray detection fraction. Very few YSOs are detected in both bands, and we 
find the radio and X-ray populations among YSOs to be largely distinct. We combine these new 
results with previous simultaneous Chandra and VLA observations of star-forming regions and 
find that YSOs with detections in both bands appear to be offset toward higher radio luminosities 
for given X-ray luminosities when compared to the Giidel-Benz relation, although even in this 
sensitive dataset most sources are too weak for the radio detections to provide information on the 
emission processes. The considerably improved sensitivity of the Expanded Very Large Array 
will provide a better census of the YSO radio population as well as better constraints on the 
emission mechanisms. 

Subject headings: stars: protostars - stars: pre-main sequence - radio continuum: stars - X-rays: stars 



1. Introduction 

Young Stellar Objects (YSOs) have been known 
for some time to be observable in X-ray and 
radio emission, tracers of high-energy processes 
(e.g., iFeigelson fc Montmerl c 19991. Briefly, ther- 
mal X-ray emission is thought to be produced by 
hot plasma heated in coronal-type activity while 
non-thermal radio emission traces magnetic fields 
that may be cospatial to and confining the X- 
ray emitting plasma. YSOs are typically grouped 
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into classes, starting with class O/I protostars as 
the earliest stages, and progressing through disk- 
dominated class II objects to basically diskless 
class III sources, also known as weak-line T Tauri 
stars. These objects are progressively less embed- 
ded, and the importance of first the protostellar 
envel ope and then the protostellar disks dimin- 



ishes (|Ladalll987t lAndre et al]|l993[) 



There are several effects which complicate a di- 
rect comparison of the high-energy processes of 
YSOs and solar-like stars. Circumstellar disks, 
a defining characteristic of protostars, may af- 
fect the overall shape of the magnetic field, for 
example by providing support for large magnetic 
structures connecting the central object and the 
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disk. Radio emission may be dominated by ei- 
ther nonthermal emission originating in magnetic 
fields, or thermal emission from the disk or stellar 
wind. The latter is usually interpreted as thermal 
bremsstrahlung due to ionized material while the 
former mos t frequently is (gyr o)synchrotron radi- 
ation (e.g., IDuI 



recent summa ry of the X-ray-radiq correlation in 
cool stars, see iForbrich et al. ( 2011 ). 



iquentiy is ( gyr o ) syncnrotron raai- 
M 1 198,4 lAndrlll996l : lGiidelll20oi . 



Note that thermal emission can easily conceal non- 
thermal er nission f rom u nderneath due to its opti- 
cal depth. Andre ( 1987t ) estimate that a spherical 



ionized stellar wind of only a few 10~ Mq yr~ is 
already opaque. While it may thus seem unlikely 
that there is a chance to detect any nonthermal 
emission from protostars at all, such cases have 
been reported ( Feigelson et al. 19981: Shirley et al 
20071) . 



There are a few different methods to identify 
nonthermal radio sources in contrast to thermal 
sources. The two primary indications are polariza- 
tion and negatiye spectral indices, as deriyed from 
the two simultaneously measured radio bands. 
Here, the spectral index a is defined as oc v", 
and a negatiye spectral index a < —0.1 is an indi- 
cation ofoptically thin nonthermal emission (e.g., 
Andrei Il996l) . While a two-frequency spectral in- 
dex is not a yery reliable indicator, the detec- 
tion of polarized emission is the most unambigu- 
ous sign of non-thermal emission as obseryed with 
the VLA. Howe yer, due to the small degree of po- 
lariza tion fe.g.. lWhite et al . 1992: Feigelson et al 



1998I ). high S/N ratios are required to meaning- 
fully detect or constrain such polarization. A third 
indication is the existence of short-term yariabil- 
ity. 

While it is not obyious that thermal hot-plasma 
X-ray emission and nonthermal radio emission 
should haye an obseryable underlying connection, 
certain types of solar flares and a yariety of ac- 
tiye stars haye been shown to haye correlated X- 
ray and (nonthermal) radio luminosities accord- 
ing to the so-called Giidel-Benz (GB) relation, 
ix/iR « 10 ^^=^^ Hz, extending oyer 10 orders 
of magnitude (iGiidel fc Benj|l993t iBenz fc Giidd 



19941 : lGudell2002[ ). The explanation which these 
authors suggested for the nearly linear correlation 
between X-ray and radio luminosities relied on a 
common energy reseryoir for both plasma heat- 
ing and particle acceleration, and a similar parti- 
tioning of energy into these two processes oyer 10 
orders of magnitude in stellar luminosity. For a 



The early studies on correlated X-ray and ra- 
dio luminosities of different object classes included 
some weak-line T Tauri stars (class HI sources) as 
powerful radio and X-ray sources, eyen though the 
underlying obseryations were non-simultaneous. 
It is therefore interesting to find out whether YSOs 
in earlier eyolutionary stages, particularly class I 
protostars - strong X-ray and radio sources in 
their own right - are compatible with the GB 
relation. Since YSOs are highly yariable, this 
question is preferably addressed by simultaneous 
X-ray and radio obseryations. Until now, there 
haye been few such obseryations of YSOs. V773 
Tau, a T Tauri star, became the first YSO ob- 
seryed in this way, and radio yariability exc eeding 



the X-ray yariability was f o und (.Feigelson et al 



1994t iGuenther et aLllioOol ). iBower et all (|2003h 



reported the serendipitous simultaneous radio and 
X-ray obseryation of a spectacular radio fiare to- 
ward a weak-line T Tauri star in Orion (GMR-A). 
In the first study targeting a subset of an entire 
star-forming region, Ga gne et al. (2004) obseryed 
the p Oph cloud complex with Chandra and the 
NRAO Very Large Array. Seyeral T Tauri stars 
were detected in both wavelength ranges, but no 
clear X- ray-radio corre l ation was found. Subse- 



quently, iForbrich et al.l ([20071) succeeded in de- 
tecting several class I protostars in the CrA star- 
forming region simultaneously in the X-ray and ra- 
dio wavelength ranges. Among the YSOs in CrA 
detected in both bands, most are c ompatible with 



the G B relation. Most recently, lOsten fc Wolk 
(|200 9h report the detections in both bands of sev- 
eral T Tauri stars in the LkHa 101 cluster. Here, 
no clear X-ray-radio correlation was found. 

In an attempt to observe more protostars si- 
multaneously in the X-ray and radio wavelength 
ranges, we have observed selected regions in two 
more young clusters, IC 348 and NGC 1333. Pre- 
viously published observations carried out with 
the Spitzer Space Telescope serve as an excellent 
road map to locate YSOs in early evolutionary 
stages in these clusters. Two-frequency radio ob- 
servations help us in potentially discriminating 
between nonthermal and thermal radio sources, 
based on the value of the radio spectral indices. 
We will discuss the results in the context of previ- 
ous simultaneous radio and X-ray observations of 
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YSOs. 

1.1. The target regions: IC 348 and NGC 
1333 

Located in the Perseus Molecular Cloud, both 
IC 348 and NGC 1333 are prominent near by star- 
formi ng re gions (for recent r eview s, see iHerbstI 



2008 and 



Walawender et alj 120081 respectively, 



and iBallv et al.l 120081 ) . The most accurate dis- 
tance has been determined toward NGC 1333 by 
measuri ng the parallax of th e bright maser source 



SVS 13 (iHirotaet alJl2008h . yielding 235±18 pc. 



IC 348 is slightly mo re distant at about 260 pc 
( Lombardi et al. 2010f) . The star formation activ- 
ity in both regions has been surveyed extensively. 
Most relevant for our purposes are studies using 
the Spitzer Space Telescope to identify and clas- 
sify YSOs in early evolu tionary stages by their 



Lada et al.l 



.g es by 



infrared excess emission, 
lyzed ~300 previously known members of IC 348 
and classified them by their disk properties, dis- 
tinguishing stellar (photospheric) , anemic, and 
thick-disk spectral energy distributions (SEDs). 
The main basis for this classification is the dered- 
dened Spitzer-IRAC slope. In a follow-up paper, 
Muench et al.l (j2007() identified and analyzed new 
cluster members, including embedded class I pro- 
tostars. Their clas s II categor y corresponds to 



the thick disks from iLada et al.l . even though the 



spectral indices o f some thick disk sources from 
Lada et al. (|2Q06|) qualify as cl ass I sources ac- 



cording to I Muench et al.l pOOTl ). To use a con- 



sistent cl assification, we follow the spectral index 



limits in I Muench et al.l (|2007l ) without separat- 
ing out flat-spectrum s ources . The class III ob- 



jects in iMuench et al.l (|2007f ) corresp ond to the 



anemi c and photospheric sources in iLada et al 
(|2006l ). The NGC 1333 reg i on wa s similarly sur 



veyed bv iCutermuth et al-l l 20081 ) who employ a 
three-phase technique to identify YSOs. Our focus 
here is on the youngest evolutionary stages. This 
is due both to scientific interest and the fact that 
in both regions previous surveys have been most 
complete for class I and class II sources. However, 
our study also reveals class III sources. Since the 
Spitzer studies cannot distinguish class and class 
I s ources, we ad d itiona lly use the sources identified 
bv lEnoch et al. ( 20091 ) in Perseus to find the class 
sources among the Spitzer class I sources. Com- 
pared to star-forming regions that have previously 



been targets of simultaneous X-ray and radio ob- 
servations (see above), IC 348 and NGC 1333 are 
about twice as distant as Taurus, p Oph, and CrA, 
but at less than half the distance of LkHa 101. 

1.2. Previous radio observations 

Both IC 348 and NGC 1333 have been observed 
with the NRACQ Very Large Array (VLA) prior 
to our own observations. Our regi on of inter- 
est in IC 348 has been observed by lAvila et"al 



(|200lt ) who present a deep radio dataset centered 
on HH211. Their concatenated dataset has a 
quoted rms noise sensitivity of 12 ^Jy and con- 
sists of a short C array observation obtained in 
1994 and 9.7 h of on-source time using 11 VLA 
antennas in 1999. The phase center of these ob- 
servations lies ~ 9' off our position 1 (see below) 
and ~ 2' off our position 2. 

NGC 1333 was studied in VLA observations 
of SVS 13 and its surroundings. After an ini- 
tial study (jRodriguez et al.l ll997l ) with two ob- 
servations reaching sensitivities of 30 /xJy and 
20 iiJy, deep concatenated archival datasets in the 
3.6 cm (X) and 6 cm ban ds (C) were presented 
by ( Rodriguez et al. 1999f ). While these datasets 
contain data obtained in various configurations, 
a common subset of the uv plane has been used 
for the combined data, resulting in angular reso- 
lutions of 4"~5" in both bands. The phase centers 
of these observations coincide with our southern 
pointing (see below) whereas our northern point- 
ing is off by ^ 7'. 

2. Observations and data analysis 

Both IC 348 and NGC 1333 have been observed 
with very similar observational setups. Since the 
Chandra-AGIS field of view (17' x 17') is consider- 
ably larger than the VLA primary beam, we have 
selected two VLA pointings within the Chandra 
fields (see Figure [T|). Using the VLA with two dis- 
tinct subarrays allows us to record two different 
radio frequencies simultaneously. 



^The National Radio Astronomy Observatory is a facility of 
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2.1. IC 348 

A region centered on RA 03h43ni59.90s, Dec 
+31°58'21.70" was observed with Chandra ACIS- 
I in two 40 ksec observations on March 13 and 
on March 18, 2008 (IDs 8933 and 8944). Si- 
multaneous radio observations were carried out 
with the VLA in its C configuration (project ID 
S9056). Two subarrays were used to cover the X 
and C bands simultaneously, using 13 antennas in 
X band and 14 antennas in C band. Two point- 
ing positions within the ACIS-I field of view were 
observed alternatingly, RA 03h44ml3.58s, Dec 
-K32°00'24.21" (position 1) and RA 03h43m52.96s, 
Dec -h32°02'44.21" (position 2). Total on-source 
time for each pointing and frequency is approxi- 
mately 10.0 hours. The full- width half-maximum 
primary beam sizes in the X and C bands are 5.3' 
and 9.6', respectively. An absolute fiux density 
scale was established by an observation of 3C147 
on the second day of observations. The phase cal- 
ibrator 0336-1-323 was observed after each 20 min 
scan on one of the two pointing positions. The 
primary calibrator (3C147) was at 4.68 Jy (X- 
band) and 7.86 Jy (C-band), and the phase cal- 
ibrator 0336-1-323 had bootstrapped flux densities 
of 1.24 Jy (X-band) and 1.68 Jy (C-band). 

2.2. NGC 1333 

A region centered on RA 03h29m02.00s, Dec 
-H31°20'54.00" was observed with Chandra ACIS- 
I in two 40 ksec observations on July 5 and on July 
11, 2006 (IDs 6436 and 6437). Simultaneous ra- 
dio observations were carried out with the NRAO 
Very Large Array (VLA) in its B configuration 
(project ID S7874). Two subarrays with 13 an- 
tennas each were used to cover the X and C bands 
simultaneously. Two pointing positions within the 
ACIS-I field of view were observed alternatingly, 
RA 03h29ml4.072s, Dec -H31°22'40.06" (position 
N) and RA 03h29m04.941s, Dec -|-31°15'54.37" 
(position S). The total on-source time for each 
pointing and frequency is approximately 9.4 
hours and the typical on-source scan-duration was 
15 minutes. An absolute flux density scale was es- 
tablished by observations of 3C286 on both days 
of observations. The phase calibrator 0336-f323, 
the same that was also used for IC 348, was ob- 
served after each scan on one of the two pointing 
positions. The primary calibrator (3C286) was 



at 5.20 Jy (X-band) and 7.49 (C-band), and the 
phase calibrator 0336-1-323 had bootstrapped flux 
densities of 1.91 Jy (X-band) and 2.28 Jy (C- 
band) . 

2.3. Sensitivity 

Given the different sizes of the fields of view 
of Chandra- AGIS and the VLA, there is a trade- 
off between increasing spatial coverage by using 
several VLA pointings and increasing sensitivity 
by observing only one position. Additionally, ob- 
servations at only one frequency would be more 
sensitive than the use of two frequencies, either by 
using two subarrays or by observing separately. 
As a compromise, we opted for a setup using two 
pointings and two frequency bands: a loss of sen- 
sitivity. Compared to a setup where we would 
maximize sensitivity by observing a single posi- 
tion with all antennas at one frequency, we lose 
a factor of 2--\/2=2.8 in flux density sensitivity. 
Primarily due to the different distances to star- 
forming regions, this study is not as sen sitive as 
the p Oph st udy by Gaene et al. ( 20041 ) or the 
CrA study bv iForbrich et all (|2007t ). Our obser- 
vation strate gy is a comp r omise between the one 
employed by ICagne et al.l (|2004l ) who used a mo- 
saic of_\Tj^^_2ointings_to cover the Chandra held 
and iForbrich et al. 1 (l2006l) who maximized sensi- 
tivity by using only one pointing position and fre- 
quency. We show a more detailed comparison of 
the various simultaneous Chandra and VLA ob- 
servations of star- forming regions in Table [1] 

2.4. Data reduction and analysis 

All VLA observations were analyzed using the 
NRAO Astronomical Image Processing System 
(AIPS), following standard procedure. After data 
flagging and calibration, the primary beam ar- 
eas of both experiments were imaged using a ro- 
bustness power of for uv weighting as a com- 
promise between sensitivity and synthesized beam 
size. For the best sensitivity in subsequent source 
detection both observing epochs of the two clus- 
ters were combined for imaging. Source detection 
was performed using the AIPS task 'SAD' above a 
S/N ratio of 5. All maps have rms noise levels near 
the theoretical expectation of 18 /iJy, based on the 
VLA Exposure CalculatoiH, assuming 10 hours of 



^http : //www . via . nrao . edu/astro/gTildes/exposure/calc . html | 
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Table 1 

Sensitivities of simultaneous Chandra and VLA observations of star-forming regions 



Target 


Date 


d 
(pc) 


(ergs^^ Hz^^) 


Lx{5 cuts) 
(ergs-i) 


References 


p Oph 


May 2000 


130^ 


3.2 X 10^5 


7.8 X 10^6 


GaEue et al. (2004) 


LkHalOl 


Mar 2005 


510 


1.5 X 10^6 


1.5 X 10^8 


Osten & Wolk f2009) 


CrA 


Aug 2005 


130 


1.7 X 10^5 


1.0 X 10^7 


Forbrich et al. (2007)'' 


NGC 1333 


Jul 2006 


235 


6.9 X 10^5 


3.2 X 10^7 


this work 


IC 348 


Mar 2008 


250 


5.7 X 10^5 


3.6 X 10^7 


this work 



Note. — The radio and X-ray sensitivities have been determined using the 5a radio detection limit in 
either C or X band and assuming a 5 count detection limit of an unabsorbed source with an APE C hot- 
plasm a spectrum (with a temperature of 1 keV) in the Chandra data, as determined using PIMMS (jMukai 
19931 ). 



Gagne et al.l (j2004f ) assume d = 165 pc, but see discussion in lWilking et al.l (j200 



^From simultaneous obs ervations. For deeper radio and X-ray data, see iForbrich et al.l (|2006l ) and 
Forbrich fc PreibischI (|2007t ). 



on-source time at full bandwidth, effectively 86 
MHz, and 13 antennas. The fitted source flux 
densities were corrected for primary beam attenu- 
ation and the peak flux densities from the fits are 
reported in the tables. Radio upper limits were de- 
termined on maps that were previously corrected 
for primary-beam attenuation. The synthesized 
beam sizes are on the order of 2.5" (IC 348) and 
0.8" (NGC 1333) in X band and on the order of 
4.0" (IC 348) and 1.2" (NGC 1333) in C band. 

The Chandra X-ray observations have been re- 
duced using the ANCHORS pipeline (AN archive 
of Chandra Observations of Regions of Star forma- 
tiord). The X-ray observations of NGC 1333 have 



alrea dy been published elsewhere ([Winston et al 
I2OIOI) . Full details of the ANCHORS processing 
are also given there. For previous X-ray r esults on 
NGC 1333, see also lCetman et al.l (|2002l ). 



3http://cxc.harvard.edu/ANCH0RS 



3. Results 

3.1. Radio data 

3.1.1. IC 348 

Only five radio sources were detected in the two 
integrated IC 348 pointings above a significance 
of 5(7 (map rms 14 /iJy). In addition to using the 
AIPS task SAD, wc have manu ally searched about 
160 YSO positions repor ted by Lada et al. ( 20061 ) 
and lMuench et al.l (|20Q7t ) for lower-significance de- 
tections in either the X or C-bands, if they were 
reasonably close to the pointing centers. This 
search resulted in two additional detections, VLA 
2 and LRL 51. While VLA 2 has been reported as 
a radio source before, LRL 51 has not. Since LRL 
51 is detected in both radio bands, it seems un- 
likely that this would be due to a chance alignment 
with a noise peak; see the discussion of chance 
alignments for NGC 1333 below. 

In Table |2l we list the total of seven radio 
detections in the nominal primary beam areas 
(FWHM) from both detection methods. Five 
sources ca n be identified wit h so urces from the 
sampl es of iLada et al. and iMuench et al 

(|2007t) . These are HH211-mm, a class proto- 
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star, LRL 52590 and LRL 51, class I sources, LRL 
13 (class II), as well as LRL 49, a source listed as 
having a stellar SED (class III). Additionally, we 
note weak outflow-like emission on opposing sides 
of LRL 54459/54460 and LRL 57025, but these 
tentative detections are not listed in the source 
table. Source 7 does not have infrared or optical 
counterparts; it is thus probably extragalactic. 

To identifiy signs of nonthermal emission, we 
have also searched Stokes- V maps for detections of 
circularly polarized soources, a telltale sign for gy- 
rosynchrotron radiation. No sources were detected 
above 5a in a Stokes- map; the 3ct upper limits 
for the two brightest YSOs are 45 /iJy, correspond- 
ing to upper limits for the circular polarization of 
LRL 52590 and LRL 49 of 3.5% and 8.2%. Ad- 
ditionally, we are particularly interested in radio 
counterparts with negative spectral indices, again 
indicative of nonthermal emission. However, we 
do not find clear cases of negative spectral indices. 
While LRL 13 only has a C band counterpart and 
remains undetected in the X band, the detection is 
too weak for the upper limit of the spectral index 
to be significant. 

The low number of detections correspond to low 
detection rates since considerably more YSOs are 
located in the primary beam areas. Table |3] lists 
the different source types located in the primary 
beam areas of the X-band and C-band observa- 
tions as well as the respective radio detections. As 
noted above, the C-band primary beam is larger 
than the X-band primary beam. All YSO classes 
have very low detection rates. Overall, only 10% 
of the class O-II YSOs in the C-band primary beam 
areas were detected. 

Our list of detec tions contains two of the four 
sources reported bv lAvila et al. (2001) in a deeper 
observation. For the two undetected sources, we 
can report upper limits. VLA 1 was reported with 
an X-band flux density of 0.14 mJy where we find 
an upper limit of < 0.07 mJy {5a), and for VLA 4, 
the reported X-band flux density of 0.59 mJy com- 
pares to our upper limit of < 0.08 mJy {5a). 
These two sources appear to have no known coun- 
terparts at other wavelengths and probably are 
variable extragalactic radio sources. Both VLA 2 
and VLA 3 also show variability on the time scales 
of years. VLA 2 is slightly weaker, but statisti- 
cally marginally compatible with the previous re- 
sult while VLA 3 has become brighter. We sum- 



marize this comparison in Table ID Among our 
radio detections, only one source is located in the 
primary beam area of the Avila et "al] ( 2001 ) ob- 
servations, i.e., source 4, or LRL 13. However, we 



only d etect the source in C band while I Avila et al 
( 200l[ ) only discuss X band data. 



When imaging the two epochs of this experi- 
ment separately, some variability on time scales of 
days becomes apparent. Only two of the sources 
seen in the full dataset are detected at > 5a in the 
single epochs and are listed in Table [5] Interest- 
ingly, we also detect and list two more sources that 
are detected above the same threshold in only one 
of the two epochs and remain below the 5a cut- 
off in the combined data (note that these are not 
included in Table |3| . One is LRL 1888, listed by 
Muench et al. (20071) as a candidate class HI mem- 



ber, and the other is unidentified, though within 
1" of a previously listed infrared source ([PSZ2003] 
J034403.6-H320520). Both detections are too weak 
for reliable determinations of spectral indices. 

3.1.2. NGC 1333 

We have used the same methods that were used 
for IC 348 for source detection in NGC 1333. An 
automatic search above a significance cutoff of 
5a was followed by a manual search for lower- 
significance det ections toward previou sly reported 
radio sources (iRqdrfguez et al. I [l999h and YSOs 
( Gutermuth et al.[ 2008[) : the typical map rms is 
20 ^Jy. Nine sources were detected by the au- 
tomatic search while ten more sources were found 
manually. The resulting source list with 19 sources 
is shown in Table IHl 

At detection significance levels of 3cr, used in 
the manual search, the number of false detections 
due to noise peaks is no longer negligible. When 
using the total number of sources that we searched 
manually and assuming Gaussian noise, we expect 
~l-3 false detections in the X and C bands due to 
chance alignments of candidate sources and posi- 
tive image pixels with > Za significance within the 
area of a synthesized beam. Note that this only 
holds for detections in just one of the two bands. 
While Table [5] also contains sources with single- 
band detections at levels of < 5a, these could thus 
be chance alignments. 

In addition to such alignments with noise peaks, 
there are also detections of extragalactic sources in 
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Table 2 
ic 348: radio detections 



no. Pos 


RA, Dec 


Sx 


Sc 


a 


ID 




(J2000.0) 


(mJy) 


(mJy) 


(X-C) 


(lit.) 



1 2 3:43:56.816 +32:00:50.06 0.087±0.022 0.040±0.016 VLA2, HH211-nim (0) 

2 2 3:43:57.099 +32:03:03.79 0.079 + 0.015 < 0.07 (Sct) 

3 2 3:43:57.610 +32:01:37.39 0.552 + 0.017 0.542 + 0.015 0.03 + 0.11 VLA3, LRL49 (STAR) 

4 2 3:43:59.707 +32:01:53.43 < 0.08 (5cr) 0.087 + 0.014 LRL 13 (THICKJI) 

5 1 3:44:12.990 +32:01:35.71 0.057 + 0.018 0.054 + 0.017 LRL 51 (THICK,!) 

6 1 3:44:20.377 +32:01:58.60 1.312 + 0.021 1.146 + 0.019 0.25 + 0.06 LRL 52590 (I) 

7 1 3:44:31.459 +32:00:39.96 out of pb 0.227 + 0.026 - no IR counterpart 



Note. — The spectral index is only given when the corresponding error is < 0.2. The second column 
refers to a detecti on in either of the two pointings, as defined in the text. The V LA designa tions are from 



Avila et al.l (|200lh and the LRL numbers and classification s are from iLada et al 



and Muench et al 



(|2007l ). The two thick disk sources from |Lada et al.l (|2006l ) have d ifferent classes accordi ng to the spectral 



index criterion used by I Muench et al. I (l2007l) . For HH 211-mm, see iFroebrich et al.l I (l2003l) . 



Table 3 

IC 348: Radio - detected members 





FOV (C) 


FOV (X) 


det. in X or C 


Class 


3 


3 


1 


Class I 


14 


14 


2 


Class IP 


22 


12 


1 


Class III'' 


28 


13 


1 



Note. — Based on Lada et al. ( 20061 ) and Muench et al. (2007), refined with class sources from 
Enoch et al.l (|2009[) . 



^including SED type THICK from lLada et al.l (|2006l ). apart from sources LRL 51 and LRL 276 which are 
class I, see text. 



'including SED types ANEMIC and STAR from lLada et al.l (|2006D 
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Table 4 

ic 348: detections of previously known radio sources 



Src 


Position 


Sx 


Sx 


HH 211- 


(J2000) 


(Avila et al. 2001) 


(this work) 


VLA 1 


03 43 47.58 +31 59 41.45 


0.14 mJy 


< 0.07 mJy {5a) 


VLA 2 


03 43 56.80 +32 00 50.41 


0.09 mJy 


0.040 + 0.016 mJy 


VLA 3 


03 43 57.60 +32 01 37.52 


0.41 mJy 


0.542 + 0.015 mJy 


VLA 4 


03 44 12.10 +31 58 29.32 


0.59 mJy 


< 0.08 mJy {5a) 



Note. — The source positions are from lAvila et al.l (|2001l ) but have been converted to J2000 



Table 5 

IC 348: RADIO detections by epoch 



no. 


Pos 


ID, RA/Dec 


Epochs 1&2 
Sx (mJy) 


Sc (mJy) 


a 


3 


2 


LRL 49 (STAR) 


Ep 1 


0.829 + 0.025 


0.726+ 0.022 


0.24 + 0.11 








Ep 2 


0.364 + 0.024 


0.417+0.019 


-0.25 + 0.20 


6 


1 


LRL 52590 (I) 


Ep 1 


1.472 + 0.033 


1.257+ 0.028 


0.29 + 0.08 








Ep 2 


1.199 + 0.030 


1.095 + 0.023 


0.17 + 0.08 




1 


3:44:19.797 +31:59:18.63 


Ep 1 


< 0.106 {5a) 


< 0.127 {5a) 








(LRL 1888, III) 


Ep 2 


0.175 + 0.025 


0.133 + 0.022 






2 


3:44:03.652 +32:05:20.67 


Ep 1 


< 0.080 {5a) 


0.179 + 0.031 










Ep 2 


< 0.097 {5a) 


< 0.080 {5a) 





Note. — Source numbers refer to Table[2] The spectral index is only given when the corresponding error 
is < 0.2. 
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the field of view of both clusters. Above the re- 
spective 5(7 detection limits, we can expect about 
one extragalactic source in each X-band primary 
beam area and about three extra galactic sources in 



each C-band primary beam area ([Windhorst et al 



19931 ) . The probability of a chance alignment of a 
known YSO in the manual search with an extra- 
galactic source is, however, exceedingly low. 

Out of 1 9 sources, ten can be id entified with 
YSOs from iGutermuth et all (|2Q08f ). One addi- 



tional source, VLA 2a/b, has a millimeter coun- 
terpar t has been descri b ed by Chandler fc Richeil 
(|2000l ) and IChen et all (|2009l ) as a class proto- 
star, but it is an unresolved binary in our data. 
The source does not appear to have a mid-infrared 
( S pitzer-IRAC) coun t erpar t and was not reported 
by IGutermuth et al.l (|2008t ). Two of these YSOs 
have clearly negative spectral indices, the un- 
resolved source VLA 2a/b and [GMM2008] 97. 
Source 10, while not selected as a Spitzer YSO, 
has been des c ribed as an M2.7 dwarf star by 
Wilking et al.l (|2004r ). Stokes- ^ maps show that 



none of the NGC 1333 sources reported here have 
significant circular polarization. 



Previouslv. lRodrfguez et al.l (jl997D reported re- 
sults from a comparable VLA observation of the 
region, resulting in the detection of sources VLA 
1-4, all of which we can confirm, although the fiux 
densities appear to have changed (Table [7]). Addi- 
tionally, these detections can be compared to the 
deep r adio study of the region by iRodriguez et al. 
()l999l) . Their map is twice as deep as ours, with 



an rms noise of 10 fi3y, and they detect 35 sources 
above a significance of 5 tr. These sources in- 
clude 16 X-band detections above 0.10 mJy, cor- 
responding to the 5cr cutoff used in our data. All 
but one of these sources are in the primary beam 
area (FWHM) of our southern po inting, but many 



remain undetected. However, iRodrfguez et al 
(|l999f ) use a concatenated dataset consisting of 
several VLA observations in different configura- 
tions, requiring the use a common subset in the 
uv plane so that the beam size is ~ 5", emphasiz- 
ing more extended structures than our setup. We 
ha ve searched all 20 X-b and detections reported 
bv IRodrfguez et aL I (|1999[ ) for counterparts in our 
X-band data, and we find five additional radio 
sources that did not make our SNR cutoff; see Ta- 
ble [7] where we also compare the flux densities t o 



the values given bv IRodrfguez et aL ( 1997 . 19991 ) 



Only one source has a fiux density that is sta- 
tistically compatible with the previously reported 
value (within la). 

Finally, we report the radio detections of the 
different YSO evolutionary classes in Table E] The 
table lists the YSOs located in the primary beam 
areas and the detection counts. The resulting 
detection fractions are again low, if marginally 
higher than in IC 348. Across classes 0~II, 13% 
of the YSOs located in the C-band primary beam 
areas have been detected in either radio band. 

We have again checked for source variability 
by separately imaging the two epochs. Only two 
sources are detected above a significance of 5a in 
the two individual epochs, and they are listed in 
TablelH Both sources, VLA 2a/b and [GMM2008] 
97 show some variability and have consistently 
negative spectral indices. 

3.1.3. Summary 

Between the two clusters IC 348 and NGC 1333, 
we report the detections of radio counterparts to- 
ward 16 YSOs, 11 of these are located in NGC 
1333. None of these could be unambiguously iden- 
tified as nonthermal radio sources due to polariza- 
tion. Using simultaneously measured negative ra- 
dio spectral indices and variability instead as prox- 
ies for nonthermal emission yields two nonthermal 
radio counterparts in NGC 1333 (VLA 2a/b and 
[GMM2008] 97). Vice versa, only a single source, 
the class I source LRL 52590, clearly appears to be 
a thermal source when judged by its radio spectral 
index of a = 0.25 ± 0.06. 

3.2. X-ray data 

The Chandra X-ray observations cover a con- 
siderably larger area than the VLA observations. 
In the following, we mainly focus on the Chandra 
X-ray data covering the VLA primary beam areas. 

3.2.1. IC 348 

In order to reduce and process the X-ray data, 
we have used the ANCHORS pipeline. Briefiy, the 
two observations are combined for the purpose of 
detecting point sources. Two iterations with the 
ciao wavelet-based tool wavdetect were performed. 
The first pass took a square region 15' on each side 
and centered at the aimpoint at full resolution. 
The second pass examined a larger (25'x25') area 
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Table 6 
NGC 1333: radio detections 



no. 


i\T /c; 








f Tonnn r\\ 

[o ZUUU.U ) 






(mJy) 




[mjy) 




TD 
[lit.) 


8 


S 


3 


28 


57 


382 +31:14 


15 


79 


< 0.11 (5cr) 





164+ 0.025 




[G08] 3 (P) 


9 


s 


3 


28 


57 


650 +31:15 


31 


33 


0.307 + 0.026 





483 + 0.023 


-0.83 + 0.24 


VLA 1, extragal. 


10 


s 


3 


29 


00 


253 +31:13 


39 


19 


0.170 + 0.037 




< 0.14 (5cr) 




VLA13, (M2.7'') 


11 


s 


3 


29 


01 


964 +31:15 


38 


04 


1.073 + 0.021 


1 


237+ 0.021 


-0.26 + 0.07 


VLA 2a/b, MMS 3 


12 


s 


3 


29 


01 


970 +31:15 


36 


22 


0.093 + 0.021 





185 + 0.022 






13 


s 


3 


29 


03 


082 +31:15 


51 


66 


0.078 + 0.020 




< 0.12 (5cr) 




VLA 17 


14 


s 


3 


29 


03 


372 +31:16 


01 


60 


0.100 + 0.021 





122 + 0.021 




VLA 3 


15 


s 


3 


29 


03 


756 +31:16 


03 


88 


0.164 + 0.020 





107+ 0.025 




VLA 4, [G08] 29 (I) 


16 


s 


3 


29 


05 


742 +31:16 


39 


71 


0.077 + 0.021 




< 0.11 (5(7) 




VLA22JG08] 84 (II 


17 


N 


3 


29 


07 


772 +31:21 


57 


03 


0.130 + 0.019 




< 0.12 (5cr) 




[G08] 31 (I) 


18 


N 


3 


29 


10 


373 +31:21 


58 


94 


0.442 + 0.020 





573 + 0.022 


-0.48 + 0.15 


[G08] 97 (II) 


19 


s 


3 


29 


10 


420 +31:13 


32 


09 


0.179 + 0.041 




< 0.13 (5(t) 




VLA 25, [G08]6(0>^ 
[G08] 7 (0'^ ) 


20 


s 


3 


29 


11 


250 +31:18 


30 


99 


< 0.26 (5cr) 





171 + 0.028 




21 


N 


3 


29 


17 


678 +31:22 


45 


49 


0.055 + 0.018 





077+ 0.022 




[G08] 111 (II) 


22 


s 


3 


29 


20 


643 +31:15 


49 


59 


out of pb 





816 + 0.029 




VLA 32, extragal. 


23 


s 


3 


29 


20 


908 +31:15 


49 


36 


out of pb 





169 + 0.029 






24 


N 


3 


29 


21 


328 +31:23 


46 


20 


< 0.14 (5cr) 





089 + 0.024 




[G08] 117 (II) 


25 


s 


3 


29 


22 


266 +31:13 


54 


40 


< 0.77 (5cr) 





178 + 0.036 




[G08] 38 (I) 


26 


N 


3 


29 


30 


923 +31:22 


11 


82 


out of pb 





479 ± 0.033 







Note. — Running number continued from Table [2] The second column refers to a detection in cither the 
northern (N) or the southern (S) pointing. The spectral index is only given when the cor responding error 
is < 0.3. [GMM2008] in ID column abbreviated to [G08 ] ; the VLA designat i ons a r e from iRodrieuez et al 



(|l997ill999l) . For the evolutionary state of source 11, see lChandler fc Richeil (|2000[) : IChen et all 



'^from lWilking et al.l (|2004D 
'^from lEnoch et all (|2009[ ). 
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Table 7 

NGC 1333: detections of previously known X-band radio sources 



Src 


Position 




S'x (total) 


Sx (peak) 


S'c (total) 


S'c (peak) 


HH 7-11- 


(J2000) 




(mJy)^ 


(mjy, this work) 


(mJy)'^ 


(mJy, this work) 


VLA 1 


03:28:57.65 - 


f31 


15:31.5 


0.52 


0.307 ± 0.026 


1.05 


0.483 ±0.023 


VLA 2a 


03:29:01.95 - 


f31 


15:38.3 


1.43 


1.073 ± 0.021 


1.36 


1.237 ±0.021 


VLA 3 


03:29:03.37 - 


f31 


16:01.9 


0.27 


0.100 ±0.022 


0.12 


0.122 ± 0.021 


VLA 4 


03:29:03.75 - 


f31 


16:04.1 


0.32 


0.164 ±0.020 


0.21 


0.107 ±0.025 


VLA 13 


03:29:00.29 - 


f31 


13:38.6 


0.13 


0.170 ±0.037 


< 0.06 


< 0.14 (5cr) 


VLA 17 


03:29:03.08 - 


f31 


15:52.5 


0.08 


0.078 ±0.020 


< 0.04 


< 0.12 (5cr) 


VLA 22 


03:29:05.74 - 


f31 


16:39.5 


0.17 


0.073 ±0.020 


0.17 


< 0.11 (5cr) 


VLA 25 


03:29:10.46 - 


f31 


13:32.1 


0.49 


0.179 ±0.041 


0.36 


< 0.13 (5cr) 


VLA 32 


03:29:20.67 +31 


15:49.6 


2.55 


out of pb 


3.64 


0.816 ±0.029 



Note. — The source positions are from Rodriguez et al. (1997 ) and Rodriguez et al.l (|l999h but have been 
converted to J2000. Note that the flux densities from jRodri guez et all "9991 ) are total fluxes while we quote 



peak flux densities. Flux densities have been corrected for primary beam attenuation. 

1999h . re 



from Rodriguez et al. (1997 



Table 8 

NGC 1333: Radio - detected members 





FOV (C) 


FOV (X) 


det. 


Class 


9 


5 


3 


Class I 


20 


8 


4 


Class II 


57 


21 


4 


Class II/III 


2 


1 






Note. — YSOs from ICutermuth et al.l ( 2008[) . refined with class sources from lEnoch et al.l ( 2009t ) 



Table 9 

NGC 1333: radio detections by epoch 



no. Pos 



ID 



Sx (mJy) 



Sc (mJy) 



a 



11 S VLA 2a/b 

18 N [GMM2008] 97 (II) 



Ep 1 
Ep 2 
Ep 1 
Ep 2 



1.216 ±0.034 1.228 ± 0.030 -0.02 ±0.10 

1.021 ±0.026 1.236 ± 0.026 -0.35 ± 0.09 

0.613 ±0.034 0.699 ± 0.030 -0.24 ±0.18 

0.336 ±0.024 0.450 ± 0.026 -0.54 ±0.24 
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including all six ACIS chips with pixels binned by 
two. To measure source parameters, the two ob- 
servations were analyzed separately. For the ex- 
traction regions of each source, the parameters of 
the ellipse-shaped point-spre ad function are in ter- 
polated from a lookup table ([Allen et al.ll2004l ) for 
a 95% encircled energy radius. The total counts 
are then scaled to 100%. 

In the combined dataset, 149 sources were iden- 
tified. In the two single datasets, 134 sources and 
109 X-ray sources were found in datasets 8933 
and 8944, respectively. The field of view of the 
Chandra o b servat ions covers 60 YSOs listed in 
Lada et all (l2006l Tab le 2) and also 30 YSOs 
fisted in lMuench et al.1 (|20Q7l Tables 1, 2, and 11). 



Out of this total of 90 YSOs, including 5 candidate 
class III sources, 33 sources have X-ray counter- 
parts within a search radius of 2". Table [TU] lists 
the statistics of X-ray detections as a functions of 
object types, concatenated into totals for the dif- 
ferent YSO classes. The table also contains infor- 
mation on the detection statistics for the nominal 
primary beam areas of the VLA observations. 

The X-ray detection fraction for the YSOs is 
~l/3 in object classes I-III. None of the class 
sources were detected in X-rays. Table [TT] lists the 
X-ray sources with counterparts in the YSO sam- 
ples. Listed are the positional difference of the X- 
ray detections compared to the sample, IDs, X-ray 
counts in both observations and information on 
variability, spectral and SED types. Apart from 
the 33 above-mentioned sources, the table lists 



two ad ditional stellar sources listed bv lLada et al 
([2006|) without an evolutionary class. LRL 1916 



a class I source, shows the most spectacular vari- 
ability between the two observations. The source 
is detected at 75 counts in the first observation, 
then drops to a nominal single count in the sec- 
ond observation. 

3.2.2. NGC 1333 

The data reduct ion of the NGC 133 3 X-ray 
data is described bv lWinston et al. ( 2010[ ). These 
authors find X-ray counterparts to 54 previously 
known YSOs and identify 41 additional X-ray 
sources with infrared counterparts as diskless 
class III YSOs. The field of view of the Chan- 
dra observations cov ers 123 out o f 137 candidate 
YSOs list ed bv .Gutermuth et al.l IIOOS). As dis- 
cussed by Winston et al.l (|201Clh . roughly a third 



of the class I and flat-spectrum protostars and 
half of the class II sources are detected in X-rays. 
Again, none of the class sources have X-ray 
counterparts. 

4. Discussion 

4.1. Combined X-ray and radio properties 

The combination of the X-ray and radio data 
has been performed with a search radius of 1". 
Using larger search radii up to 4" did not increase 
the number of sources detected in both bands. 
In total, eight YSOs have been detected in both 
bands, four in each cluster. The X-ray fit results 
from ANCHORS, assuming an absorbed single- 
temperature APEC plasma, as well as X-ray and 
radio luminosities for these sources are shown in 
Table [TU The X-ray and radio source positions 
coincide within the errors. 

While the VLA and Chandra have similar an- 
gular resolution, the existence of multiple sources 
within a scale of about 1" cannot be constrained 
from such data. Multiplicity could be due to both 
actual multiple stellar systems and a single ob- 
ject with distinct jet components. As a result, a 
perceived mismatch of X-ray and radio properties 
could be due to an unresolved multiple source. 

Using the X-ray spectral parameters to predict 
the level of optically thin thermal free-free radio 
emission returns predicted flux densities that are 
up to 10,000 times less than what is observed. 
Other than these sources, several more are de- 
tected in just one of the two bands, but the over- 
whelming majority of Spzteer-selected YSOs are 
neither detected in X-rays, nor in centimetric ra- 
dio emission. 

In the following, we will separately assess the 
radio and X-ray detection fractions for class I and 
class II YSOs in the two clusters. For the sake 
of this discussion, we will leave out the class III 
sources since they are less completely catalogued 
in both clusters. In order to obtain a comparison 
between these detection fractions, we only take 
into account sources in the VLA primary beam 
(FWHM) areas. As shown in Table [Til separately 
for IC348 and NGC 1333, the detection fractions 
are very diverse. Note that some of the fractions 
are based on very low numbers and thus are quite 
uncertain. The radio detection fractions decline 
with evolutionary stage, similarly in both clusters. 
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Table 10 
ic 348: x-ray - detected ysos 





in FOV-CXO 


w/X-ray 


also in FOV-VLA(C) 


w/X-ray 


Class 


3 





3 





Class I 


15 


7 


14 


5 


Class II 


29 


10 


21 


5 


Class III 


37 


16 


24 


10 



but no such trend is apparent in the X-ray detec- 
tion fractions. The X-ray detection fractions of 
class I and II YSOs are up to ten times higher than 
the radio detection fractions, but there is also the 
curious case of class I protostars in NGC 1333 for 
which the radio and X-ray detection fractions are 
similar, but nevertheless due to almost completely 
distinct sets of sources. Only a minor fraction of 
YSOs is detected in both bands. The present re- 
sults thus do not require a tight connection be- 
tween the observed X-ray and radio emission. 

Compared to the range expected by the GB re- 
lation, the X-ray data are much more sensitive 
than the radio data, placing any radio detection 
with an X-ray upper limit off the Giidel-Benz re- 
lation by several orders of magnitude. For these 
data points, it is important to keep in mind that 
the observations in the two bands were carried 
out simultaneously so that variability cannot ex- 
plain this discrepancy. While some of the YSOs 
with radio counterparts and X-ray upper limits 
are close to the radio detection limit and may be 
chance alignments (see above) , there are also some 
clearly detected radio sources among the YSOs 
that do not have X-ray counterparts. The most 
striking example is the class protostar source 11 
in NGC 1333 with a radio flux density correspond- 
ing to a luminosity of — 7.0 x 10^^ ergs~^ but 
no X-ray detection. While here, any X-ray emis- 
sion is likely extincted, the fact that this source 
is entirely incompatible with the GB relation may 
also point to a different radio emission mechanism, 
even though source 11 is among our candidate 
nonthermal sources due to its negative spectral in- 
dex. 

Since this comparison is limited by the sensi- 
tivity of the radio data presented here, it is too 
early to draw conclusions concerning the radio 



populations of the two clusters. However, we can 
compare our result to previously published ob- 
servations. Applying the 5cr sensitivity limit of 
the present study to the results on CrA, as re- 
ported in the deep ra dio observations reported in 
Forbrich et al.l (|2006l ) and scaling to the distance 
of Perseus, shows that about half of the radio 
sources there would still be detected. While it 
thus seems likely that a genuine lack of radio ac- 
tivity plays a role in the low detection rate, the 
numbers involved are too low to warrant a firm 
conclusion. 

Leaving aside the issue of non-detections in 
one of the two bands, we can compare the YSOs 
that were detected in both bands with correspond- 
ing sources from previous experiments. In spite 
of the fact that we do not know the dominat- 
ing emission mechanism for most sources, it is 
also interesting to compare these YSOs to the 
empirical GB relation. Figure [2] shows all such 
sources fro m IC 348, NGC 1333, p Oph. CrA, and 
LkHo; 101 dCagne et al.ll2004l:lForbrich et al.ll2007 : 



Osten fc Wolkll2009[ ). again relative to the GB re- 
lation. Also, we conservatively exclude sources 
with radio detection significance levels of < Str 
or less than 100 X-ray counts. Note that since 
the radio luminosities do not depend on a spectral 
fit, they are dominated by the measurement er- 
ror. The X-ray luminosities, in contrast, can have 
uncertainties of order unity. Two things become 
apparent in this plot: 1) Compared to the GB 
relation, the YSOs appear to have higher radio 
luminosities for a given X-ray luminosity, and 2) 
there are no clear trends with evolutionary stage. 
About half of the sources are compatible with the 
GB relation even though there are sources that lie 
considerably below the GB relation, outside of its 
approximate bounds of plus or minus one order of 
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Table 11 
ic 348: x-ray-detected ysos 



LRL 


RA 


Dec 


Offset 


X-ray cnts. 


X-ray cnts. 


cnt. 


SpT 


SED type 




(Chandra) 


(J2000.0) 


(") 


(8933,raw) 


(8944,raw) 


ratio 




(IRAC) 


1840 


03:43:19.91 


+32:02:43.13 


1.75 


33 


24 


1.4 




III 


273 


03:43:52.01 


+32:03:40.45 


1.11 


20 


16 


1.3 




III 


124 


03:43:54.65 


1 or* r\r\ <a/a /^/^ 

+32:00:29.96 


0.24 


27 


42 


1.6 


M4.25 


STAR 


26 


03:43:56.06 


+32:02:13.35 


0.37 


19 


12 


1.6 


K7 


THICK 


49 


03:43:57.62 


1 o (A r\ 1 o 

+32:01:37.47 


0.33 


56 


57 


1.0 


MO. 5 


STAR 


13 


03:43:59.67 


+32:01:54.10 


0.38 


182 


107 


1.7 


MO. 5 


IHICK 


160 


03:44:02.62 


1 o o r\ 1 or" rA/^ 

+32:01:35.06 


0.37 


29 


21 


1.4 


M4.75 


STAR 


54460 


03:44:02.66 


+32:01:59.35 


0.54 


4 


11 


2.8 




I 


22021 


03:44:03.64 


1 o r\ /~\c\ oo CAP' 

+32:02:33.25 


0.31 


23 


17 


1.4 


M5 




1916 


03:44:05.79 


1 O (A /A/~v (AO P' P' 

+32:00:28.55 


0.14 


75 


1 


75.0 




I 


214 


r\o r\T r^i 

03:44:07.54 


+32:04:08.65 


0.39 


10 


17 


1.7 


M4.75 


IHICK 


276 


r\0 r\A OO 

03:44:09.22 


1 OO /AO OO 1 /i 

+32:02:38.14 


0.36 


15 


13 


1.2 


MO 


IHICK 


173 


03:44:10.23 


1 OO r\ A /AP" P'O 

+32:04:05.53 


1.62 


6 


18 


3.0 


M5.75 


THICK 


105 


03:44:11.30 


1 OO /A/^ "1 O 01 

+32:06:12.31 


0.57 


112 


83 


1.3 


MO 


STAR 


51 


03:44:13.00 


1 OO 01 OP" P* 1 

+32:01:35.51 


0.43 


14 


31 


2.2 




THICK^ 


54519 


/A *A A A -f O '~7CA 

03:44:13.78 


1 01 P'P' 0,1 OP* 

+31:55:34.85 


0.47 


92 


86 


1.1 




HI 


31 


03:44:18.21 


1 OO r\ A p" T o/^ 

+32:04:57.26 


0.65 


219 


136 


1.6 


Gl 


AMKMIC 


99 


03:44:19.33 


1 ooriTO/1 11 

+32:07:34.11 


1.12 


105 


171 


1.6 


M3.75 


IHICK 


1888 


03:44:19.78 


1 01 P'/A 10 ZAP' 

+31:59:19.05 


0.05 


46 


159 


3.5 




HI 


210 


y~i O <1 A (A/A r\ 1 

03:44:20.01 


1 OO /A/1 /IP" <I/A 

+32:06:45.49 


0.15 


45 


47 


1.0 


M3.5 


STAR 


52590 


/AO A A r\/A 

03:44:20.43 


1 OO r\i r" PT OO 

+32:01:57.93 


0.78 


50 


75 


1.5 




I 


119 


03:44:21.29 


+32:05:02.66 


0.49 


91 


153 


1.7 


M2.5 


STAR 


1889 


03:44:21.35 


+31:59:32.80 


0.10 


15 


13 


1.2 




I 


125 


03:44:21.62 


+32:06:25.71 


1.07 


70 


171 


2.4 


M2.75 


STAR 


61 


03:44:22.28 


+32:05:43.23 


0.54 


54 


80 


1.5 


K8 


THICK 


72 


03:44:22.58 


+32:01:53.75 


0.07 


81 


130 


1.6 


M2.5 


ANEMIC 


30190 


03:44:23.65 


+32:06:47.67 


0.87 


158 


158 


1.0 


M2.5 




123 


03:44:24.60 


+32:03:56.97 


0.45 


119 


231 


1.9 


Ml 


STAR 


97 


03:44:25.58 


+32:06:17.14 


0.33 


39 


699 


17.9 


M2.25 


THICK 


5 


03:44:26.05 


+32:04:30.18 


0.37 


756 


527 


1.4 


G8 


THICK 


62 


03:44:26.68 


+32:03:58.52 


0.70 


103 


96 


1.1 


M4.75 


STAR 


69 


03:44:27.05 


+32:04:43.58 


0.43 


620 


471 


1.3 


Ml 


STAR 


68 


03:44:28.49 


+31:59:54.52 


0.49 


10 


18 


1.8 


M3.5 


THICK 


55 


03:44:31.40 


+32:00:14.16 


0.79 


60 


153 


2.6 


MO. 5 


THICK 


52648 


03:44:34.44 


+31:58:00.46 


1.05 


50 


41 


1.2 




I 



No te. — Spectral types and classification information are from Lada et al. ( 2006h and Muench et al 
()2007l) . 



^class I, see text 
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Table 12 

NGC 1333: X-ray - detected class O-II YSOs in NGC 1333 





in FOV-CXO 


w/X-ray 


also in FOV-VLA(C) 


w/X-ray 


Class 


7 





7 





Class I 


26 


6 


21 


5 


Class II 


83 


44 


74 


38 



Note. — X-ray data from IWinston et al. I (I2OIOI) 



Table 13 

ysos detected in both the x-ray and the radio range 



Source Class raw cnts. n(II) kT 

(total) (10^2 cm-2) (keV) 



Lx Lfl(X-band^ ) 
erg s^i erg s"^ IIz"^ 



LRL 51 


I 


45 




















8.6 


X 


1015 


[GMM2008] 31 


I 


56 




















4.3 


X 


1015 


LRL 52590 


I 


125 


1.87 


± 


0.30 


9.21 


± 


2.10 


4.4 


X 


1029 


1.1 


X 


lOi^ 


LRL 13 


II 


289 


1.42 


± 


0.25 


2.20 


± 


0.47 


8.4 


X 


1029 


6.5 


X 


1015 


[GMM2008] 84 


II 


455 


1.36 


± 


0.15 


3.40 


± 


0.56 


7.4 


X 


1029 


5.1 


X 


1015 


[GMM2008] 97 


II 


3796 


1.03 


± 


0.04 


1.97 


± 


0.07 


5.7 


X 


1030 


3.1 


X 


1016 


[GMM2008] 111 


II 


6397 


0.28 


± 


0.01 


2.46 


± 


0.08 


5.0 


X 


1030 


3.6 


X 


1015 


LRL 49 


III 


113 


1.02 


± 


0.24 


2.44 


± 


0.72 


2.9 


X 


1029 


4.5 


X 


1016 



Note. — X-ray information is from ANCHORS for an energy range of 0.3-8.0 keV, spectral fits are only 
shown for sources with more than 100 counts. 



^C band for LRL 13 



Table 14 

Combined Radio-X-ray detection statistics (IC 348+NGC 1333) 





Cluster 


FOV (C) 


X-ray 


radio** 


radio(X or C) & X-ray 


Class 


IC 348 


3 


(0%) 


1 (33%) 


(0%) 


Class I 


IC 348 


14 


5 (36%) 


2 (14%) 


2 (14%) 


Class II 


IC 348 


21 


5 (24%) 


1 (5%) 


1 (5%) 


Class 


NGC 1333 


8 


(0%) 


3 (38%) 


(0%) 


Class I 


NGC 1333 


20 


5 (25%) 


4 (20%) 


1 (5%) 


Class II 


NGC 1333 


66 


38 (58%) 


4 (6%) 


3 (5%) 



^In C-band primary beam area. 
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50.0 40.0 30.0 20.0 10.0 3:44:00.0 50.O 40.0 30.0 20.0 43:10.0 ^ 50.0 40.0 30.0 20.0 10.0 3:29:00.0 50.0 40.0 30.0 28:20.0 

Right ascension Right ascension 

Fig. 1. — Merged Chandra pointings toward IC348 (left panel) and NGC 1333 (right panel). The eastern 
and western pointings in IC 348 are positions 1 and 2, respectively, and the northern and southern pointings 
in NGC 1333 are positions N and S. The large circles indicate the primary beam sizes for the VLA pointings, 
the C band having a larger primary beam than the X band. YSOs detected in both bands (from Table [T^ 
are marked by small circles. 



magnitude. These sources are of all three evolu- 
tionary classes. Since all of these data points are 
due to simultaneous observations, variability plays 
only a minor role. Also, there are no clear trends 
concerning the type of emission and compatibility 
with the GB relation. For example, the thermal 
source LRL 52590 is the source that is least com- 
patible with the GB relation in Figure [2] How- 
ever, also note the opposite example of source 11 
above. These differences may indicate different ra- 
dio emission mechanisms. 

In order to judge the significance of this result 
for YSOs in general it is important to keep in mind 
that this plot contains only a minor fraction of 
the YSOs that have been covered in the simulta- 
neous X-ray and radio experiments. While many 
of the sources plotted appear radio bright, many 
other sources remain undetected. The detections 
are due to sources with extreme radio and X-ray 
luminosities, even though they do not represent all 
sources in this luminosity range. To illustrate this 
point, we can use the GB relation to determine 
the X-ray luminosity that corresponds to the ap- 
proximate 5(7 radio detection limit for IC 348 and 



NGC 1333, Lx ^ 6x 10^° ergs^^. This X-ray 
luminosity limit would include the X-ray bright- 
est ~ 25% G and K pre-main seq uence stars in 
Orion (jPreibisch fc FeigelsonI 120051 ). None of the 
X-ray detected YSOs in NGC 1333 and IC 348 ac- 
tually reach this luminosity. It is, however, note- 
worthy tha t six infrared-identified YSOs in NGC 
1333 (from Gutermuth et al.ll2008 ) have X-ray lu- 
minosities up to one order of magnitude below this 
number and lie within the VLA primary beam ar- 
eas, but they do not have radio counterparts. A 
similar trend can b e seen in the data presented by 
Gagne et al.l (|2004l ). There, in p Oph, 26 X-ray 



sources are more luminous than the expectation 
for the lower limit of the GB relation at the 5ct ra- 
dio detection limit. Of these X-ray sources, only 
5 have radio counterparts. Given the radio sen- 
sitivity limits of the various studies, as shown in 
Figure [21 sources that are undetected in the radio 
range could, however, well be compatible with the 
GB relation at lower luminosities. 

A shift toward higher radio luminosities for 
gi ven X-ray luminositie s had already been noted 
bv lCiidel fc Benz (1993) who found that the most 
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1e+32 



1e+31 



1e+30 




1e+29 



IRS2 
f62624.0 ^ 

+ 



pOph-J 162603.0 

XL1-J043019.14 
>KL1-J042953.98 



L1-J043001.15 
X 1333-97 



^/kL1-J043010.87 
pOph-JI 62622.4 

+ IRS7W 



pOph-J 162643.8 

X X 

pOph-JI 62630.5 +LRL 52590 
XLRL 49 



_i_ 



1e+15 



1e+16 



1e+17 



Lr (erg/s/Hz) 



1e+18 



Fig. 2. — Plot of the X-ray and radio luminosities of all YSOs that have been simultaneously detected 
in both the X-ray and radio ranges by Chandra and the VLA, excluding sources with a radio detection 
significance of < 5ct and less than 100 counts in X-rays. The symbols '-I-', 'x', and '*' indicate YSO classes I , 
II, and III, respectively. Apart from sources reported in this pape r, and sources in p Op h from lGagne et al 



(|2004[ ). sc aled to a distance of 130 pc, IRS sources are in CrA (jForbrich et al.l |2007^ . and LI stands for 



LkHalOl (|Osten fc Wolkll2009f) . The c entral hne ma rked "GB" indicates the Giidel-Benz relation with its 
upper and lower bounds, as reported bv lGiidell (|2002l ). shown as dashed lines. The vertical lines indicate the 
5(7 sensitivities of the different observations, as converted into radio luminosities at the respective distances. 
The limits are strictly for the simultaneous observations only (some regions have deeper radio and/or X-ray 
data from separate observations). The X-ray sensitivity limits are lower than the X-ray luminosity range 
shown (see Table [1]). 
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luminous sources were slightly offset toward higher 
radio luminosities, y et still within t he qu oted lim- 
its of the relation. Guedel et al suggest 
that not only a high X-ray luminosity, but rather 
the presence of a hot plasma component is condi- 
tional for the detectability of a late-type star as a 
strong, nonthermal radio source. However, since 
we also find indications for YSOs with radio but 
no X-ray counterparts, this may not directly apply 
to YSOs. They also find that while X-ray emis- 
sion saturates that does not seem to be the case 
for the non-thermal radio emission. Consequently, 
if all sources discussed here were X-ray-saturated, 
that could also explain an offset from the GB re- 
lation. However, it seems unlikely that this would 
explain the offsets shown in Figure [2] where the 
YSOs that are furthest off the GB relation have 
the lowest X-ray luminosities. 

One important caveat in these comparisons to 
the GB relation is the nature of the radio emis- 
sion for each single detection since the correlation 
is only valid for nonthermal emission. Among the 
sources reported in this paper, no source can con- 
clusively be proven to be nonthermal due to the 
relatively low signal-to-noise ratios: No polariza- 
tion was detected, and the error bars on the ra- 
dio spectral indices also do not allow us to use 
them as definitive criteria. Among the class I 
sources in the literature, CrA-IRS 5 is known to 



show non-thermal radio emission (jFeigelson et al 



19981 ) . This source is compatible with the GB rela- 
tion (jPorbrich et all l200a [20071) . Potentially, the 



outlier sources have additional or alternative ther- 
mal radio components that would shift them to- 
ward higher radio luminosities. This thermal com- 
ponent would, however, in some cases have to ac- 
count for several orders of magnitude in luminos- 
ity. A third and again completely different possi- 
bility involves coherent e mission processes which, 
while present on the Sun (|Benz et al.ll2005ll2006[ ). 
have not yet been identified in YSOs. 

While also the thermal X-ray emission of YSOs 
is known to be due to at least two different sources 
(coronal magnetic activity and accretion) , this re- 
sults primarily in different temperature compo- 
nents in the X-ray emission, generally of similar 
luminosities. For our comparison of X-ray and ra- 
dio luminosities, the interpretation of X-ray lumi- 
nosities should thus be more straightforward. 



5. Summary and Conclusions 

We present simultaneous X-ray and radio obser- 
vations of IC 348 and NGC 1333. Only three other 
clusters have previously been observed in this way. 
In summary, our main conclusions are as follows: 

• The radio detection fraction for class I and 
class H YSOs is generally much lower than 
the X-ray detection fraction, particularly for 
the class H sources. The radio detections 
are not solely a subset of the X-ray detected 
sources. Very few YSOs (~ 5%) are detected 
in both bands. The percentage of such X- 
ray-radio detections is marginally higher for 
class I protostars than for class H sources. 

• We find a few YSOs that are detected in 
only one band (X-ray or radio), even though 
they should be easily detectable in the other 
band if the GB relation is applicable to 
these sources. Extinction helps to explain X- 
ray non-detections, but whether radio non- 
detections of X-ray sources are due to insuf- 
ficient sensitivity remains unclear. 

• The currently known sample of YSOs that 
have been simultaneously detected in both 
the X-ray and radio bands (across several 
clusters) shows that the corresponding lu- 
minosities of only about half of these sources 
are within one order of magnitude of the GB 
relation while the other half is shifted toward 
higher radio luminosities for a given X-ray 
luminosity. 

• Compared to early results on the GB rela- 
tion almost 20 years ago, simultaneous X- 
ray and radio data are now available for sev- 
eral source populations, minimizing the ef- 
fect of variability. Also, the X-ray proper- 
ties of YSOs in particular are much better 
defined while there has been comparably lit- 
tle progress on their centimetric radio prop- 
erties. 

• A key uncertainty in these comparisons is 
the nature of the radio emission since these 
relatively weak radio sources cannot unam- 
biguously be identified as nonthermal radio 
sources. Also the radio luminosity function 
for these sources is basically unknown. Ra- 
dio sensitivity strongly limits the size of the 
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datasct which can be compared to the GB 
relation. With the advent of the Expanded 
Very Large Array, the number of YSOs with 
defined radio properties can be expected to 
increase significantly. With such observa- 
tions, it will be possible to look for fainter ra- 
dio sources and to determine the radio lumi- 
nosity function for entire YSO populations. 

Support for this work was provided by the 
National Aeronautics and Space Administration 
through Chandra Award Numbers GO6-7015X 
and GO8-9017X issued by the Chandra X-ray Ob- 
servatory Center, which is operated by the Smith- 
sonian Astrophysical Observatory for and on be- 
half of the National Aeronautics Space Admin- 
istration under contract NAS8-03060. S.J.W. is 
supported by NASA contract NAS8-03060 (Chan- 
dra). The NRAO VLA data are from programs 
S7874 (NGC1333) and S9056 (IC 348). 
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